home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / claein.z / claein
Text File  |  1996-03-14  |  3KB  |  133 lines
  1.  
  2.  
  3.  
  4. CCCCLLLLAAAAEEEEIIIINNNN((((3333FFFF))))                                                          CCCCLLLLAAAAEEEEIIIINNNN((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      CLAEIN - use inverse iteration to find a right or left eigenvector
  10.      corresponding to the eigenvalue W of a complex upper Hessenberg matrix H
  11.  
  12. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  13.      SUBROUTINE CLAEIN( RIGHTV, NOINIT, N, H, LDH, W, V, B, LDB, RWORK, EPS3,
  14.                         SMLNUM, INFO )
  15.  
  16.          LOGICAL        NOINIT, RIGHTV
  17.  
  18.          INTEGER        INFO, LDB, LDH, N
  19.  
  20.          REAL           EPS3, SMLNUM
  21.  
  22.          COMPLEX        W
  23.  
  24.          REAL           RWORK( * )
  25.  
  26.          COMPLEX        B( LDB, * ), H( LDH, * ), V( * )
  27.  
  28. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  29.      CLAEIN uses inverse iteration to find a right or left eigenvector
  30.      corresponding to the eigenvalue W of a complex upper Hessenberg matrix H.
  31.  
  32.  
  33. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  34.      RIGHTV   (input) LOGICAL
  35.               = .TRUE. : compute right eigenvector;
  36.               = .FALSE.: compute left eigenvector.
  37.  
  38.      NOINIT   (input) LOGICAL
  39.               = .TRUE. : no initial vector supplied in V
  40.               = .FALSE.: initial vector supplied in V.
  41.  
  42.      N       (input) INTEGER
  43.              The order of the matrix H.  N >= 0.
  44.  
  45.      H       (input) COMPLEX array, dimension (LDH,N)
  46.              The upper Hessenberg matrix H.
  47.  
  48.      LDH     (input) INTEGER
  49.              The leading dimension of the array H.  LDH >= max(1,N).
  50.  
  51.      W       (input) COMPLEX
  52.              The eigenvalue of H whose corresponding right or left eigenvector
  53.              is to be computed.
  54.  
  55.      V       (input/output) COMPLEX array, dimension (N)
  56.              On entry, if NOINIT = .FALSE., V must contain a starting vector
  57.              for inverse iteration; otherwise V need not be set.  On exit, V
  58.              contains the computed eigenvector, normalized so that the
  59.              component of largest magnitude has magnitude 1; here the
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. CCCCLLLLAAAAEEEEIIIINNNN((((3333FFFF))))                                                          CCCCLLLLAAAAEEEEIIIINNNN((((3333FFFF))))
  71.  
  72.  
  73.  
  74.              magnitude of a complex number (x,y) is taken to be |x| + |y|.
  75.  
  76.      B       (workspace) COMPLEX array, dimension (LDB,N)
  77.  
  78.      LDB     (input) INTEGER
  79.              The leading dimension of the array B.  LDB >= max(1,N).
  80.  
  81.      RWORK   (workspace) REAL array, dimension (N)
  82.  
  83.      EPS3    (input) REAL
  84.              A small machine-dependent value which is used to perturb close
  85.              eigenvalues, and to replace zero pivots.
  86.  
  87.      SMLNUM  (input) REAL
  88.              A machine-dependent value close to the underflow threshold.
  89.  
  90.      INFO    (output) INTEGER
  91.              = 0:  successful exit
  92.              = 1:  inverse iteration did not converge; V is set to the last
  93.              iterate.
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.